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LOWER BOUNDS FOR INTERIOR NODAL SETS 
OF STEKLOV EIGENFUNCTIONS 

CHRISTOPHER D. SOGGE, XING WANG, AND JIUYI ZHU 


Abstract. We study the interior nodal sets, Z\ of Steklov eigenfunctions in an n- 
dimensional relatively compact manifolds M with boundary and show that one has 

2 —n 

the lower bounds \Zx\ > cA 2 for the size of its (n — l)-dimensional Hausdorff 
measure. The proof is based on a Dong-type identity and estimates for the gradient 
of Steklov eigenfunctions, similar to those in m and m, respectively. 


1. Introduction. 

This article is concerned with lower bounds for the size of nodal sets, 

(1.1) Zx = {x € M : exix) = 0}, 

of real Steklov eigenfunctions in a smooth relatively compact manifold (M, g) of dimension 
n>2 with boundary dM. These eigenfunctions are solutions of the equation 

Q2) fAgeA = 0, inM 

[diyCx = Xex, on (9M, 

where ly is the unit outward normal on dM. 

The Steklov eigenfunctions were introduced by Steklov [T7] in 1902. They describe 
the vibration of a free membrane with uniformly distributed mass on the boundary. The 
equation 113 ) was studied by Calderon [3] as its solutions can be regarded as eigenfunc¬ 
tions of the Dirichlet to Neumann map. 

More specifically, the ex in 113 satisfy the eigenvalue problem 

Pex = Xex, 

if the Dirichlet to Neumann operator P is defined as 

Pf = d,Hf\sM. 

where for / G C°°{dM), Hf = u is the harmonic extension of / into M, i.e., the solution 
of 

J Agu{x) =0, X € M 
|u(a;) = f{x), X G dM. 

It is well known that P is a self-adjoint classical pseudodifferential operator of order 
one whose principal symbol agrees with that of the square root of minus the boundary 
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Laplacian on dM coming from the metric 
of real eigenfunctions {ex ^} such that 

Pex^=Xjex^, and 

The spectrum Xj is discrete, with 

0 = Ao < Ai < A 2 < ..., and Xj —>■ 00 . 

Recently there has been much work on the study of nodal sets of Steklov eigenfunctions. 
It has largely been focused on the size of the nodal set 

A/a = {x e dM : ex(x) = 0} 

on the boundary dM of M. Bellova and Lin [T] proved that |A/a| < CA®, if |A/a| denotes 
d—1 dimensional Hausdorff measure with here d = n — 1 denoting the dimension of dM. 
Later, Zelditch [23] improved these results and gave the optimal upper bound |A/a| < CX 
for analytic manifolds using microlocal analysis. In the smooth case, the last two authors 
showed in [5D] showed that 

(1.3) \Mx\>cX^, 

assuming that 0 is a regular value for ex. This agrees with the best known general lower 
bounds for the boundary less case (see below), but in both [23] and [20] the nonlocal 
nature of the operators defining the eigenfunctions presented an obstacle which had to 
be overcome. 

By the maximum principle, we know that the nodal sets in M must always intersect 
the boundary dM. In other words, there can be no component of the nodal set which 
is closed in M. Thus, it is natural to study the size of the nodal set in the interior, M. 
This question was also raised by Girouard and Polterovich in [9]. 

Let us briefly review the literature concerning the study of nodal sets for compact 
boundaryless Riemannian manifolds. Let ipx denote an L^-normalized eigenfunction on 
of the Laplace-Beltrami operator on such a smooth n-dimensional manifold, i.e., 

-Ag'il^x = A^V'a- 

Yau conjectured in [^ that one should have 

cA < jZxi < CX, 

if Zx denotes the nodal set of ipX: and jZxj its (n — I)-dimensional Hausdorff measure. 
In the real analytic case both the upper and lower bounds were established by Donnelly 
and Fefferman [B]. The lower bound was established in the two-dimensional case by 
Briining [^ and Yau (unpublished); however, in all other cases, the conjecture remains 
open in the smooth case. Recently there has been much work on establishing lower 
bounds in the smooth case when n > 3. Golding and Minicozzi [4] and then later the 
first author and Zelditch [18], [19] showed that 

(1.4) |^a|>cA^, 

which matches up with the lower bounds in (11.31) which were obtained later. Another 
proof of (HU) was given by Hezari and the second author in m 


Furthermore, there is an orthonormal basis 


IdM 


ex,ex,dVg = 6j. 
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The arguments in m.m and m involved establishing a Dong-type identity, similar 
to the one in [^, and then using either lower bounds for the -norms of ip\ or upper- 
bounds for its gradient. We shall use similar arguments to establish our main result 
concerning lower bounds for the (n — l)-dimensional Hausdorff measure of the interior 
nodal sets of Steklov eigenfunctions contained in the following result. 

Theorem 1.1. Let M be a smooth relatively compact n-dimensional manifold with 
smooth boundary dM. Then there is a constant c > 0 so that 

(1.5) |^a|>cA^ 

for the [n — 1)-dimensional Hausdorff measure of the nodal sets given by (|l.ll) of the 
Steklov eigenfunctions m- 

We note that this lower bound is off by a half-power versus the best known lower 
bounds, (El), for the boundaryless case. We shall explain what accounts for this dif¬ 
ference after we complete the proof of Theorem 11.11 Also, it seems clear that in the 
two-dimensional case the lower bound (El is far from optimal since the arguments of 
Briining [2] and Yau (see also E) seem to give the optimal lower bound \Z\\ > cA using 
the fact that the nodal set must intersect any CA”^ ball in M if C is large enough (see 
e.g. [9]). 

2. An interior Dong-type identity for Steklov eigenfunctions. 

As in [18] we shall want to use the Gauss-Green formula to establish a Dong-type 
identity which we can use to prove our lower bound (ED- We shall be able to do this 
since the singular set 

S\ = {x € M : e\{x) = 0 and Ve\{x) = 0} 

is of Hausdorff codimension 2 or more, i.e., dim S\ < n — 2. This is true for S'a flM since 
e\ is harmonic in M (see e.g. (TUI Chapter 4]), while one can, for instance see that the 
same is true for S\ DdM using the doubling lemma in [53]. In addition, for each A, there 
are only finitely many nodal domains (see e.g. |9]). Consequently, we may write M as 
the (essentially) disjoint union 

k\ mx 

(2.1) M = U (D+ U Z+ U Y+) U U (D- U Z" U Y"), 

i=i j=i 

where Df and DJ are the connected components oi {x £ M : eA(a;) > 0} and {x £ M : 
e\{x) < 0}, respectively, while Z^ = dD^ n M and n dM. Thus, 

k\ mx 

2=1 i = l 




dM=\jY+ U\J Y-. 

1=1 


i=l 


and 
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The boundary of in M is U . Since S\ has codimension 2 or more and dM 
is smooth, we may use the Gauss-Green formula (see e.g. Theorem 1 on p. 209 of H) 
for any / e C°°{M) to get 


[ Agfex dV = [ fAgCx dV - [ fd„ex dS + [ d^fex 
Jd+ Jot ddD+ JdDt 


dS 


= -A 


fexdS ■ 


/|VeA|d^- 


d^fex dS. 


lYi 


Here du denotes the outward Riemann derivative on and we used the equation dO]) 

to get the last equality. Rearranging, we see from above that 


( 2 . 2 ) 


A/ fexdS- d,fexdS+ Agf exdV = f\Vex\dS. 


’Y, 


+ 

k k 


'D2 


dS 


Similarly for each negative nodal domain we have 

[ Agfex dV = [ fAgex dV - [ fd^ex dS + [ d^fex 
Jd22 Jd 2 JdD- 

= -A / _ fexdS - f _ f\Vex\dS + f _ d^fex dS, 
dY^. Jz^, Jz^, 

using in the last step that on each Z'^, unlike on each Z'^, d^ex = IVcaI since ex increases 
as it crosses from D'jf. Rearranging this time leads to 


(2.3) 


A/ fexdS- d,fexdS+ Agf exdV = - f\S/ex\dS. 


'D7 


Since Ca > 0 in and Ca < 0 in , we can combine (|2.2I) and dH into 

(2.4) x[ /|eA|d^- / 9./|eA|d^+ / Agf\ex\dV=[ f\Vex\dS. 

JYf JyS Jof Jzf 


Since almost every point in Zx belongs to exaclty one Zf and one Z^ and almost 
every point in dM belongs to just one of the sets , if we sum up the identity (12.4L 
we conclude that we have the Dong-type identity 

(2.5) \[ f\ex\dS- [ d,f\ex\dS+ j A 3 /|eA|dR = 2/ /|VeA|d^. 

JdM JdM JM JZx 

Of course if / = 1 this simplifies to 

(2.6) x[ \ex\dS = 2[ \Vex\dS, 

JdM J Zx 

which is what we shall use in our proof of Theorem ll.il 


3. Interior estimates for Steklov eigenfunctions. 

We shall prove interior estimates for the Ca which are natural analogs of the ones 
obtained earlier in the boundaryless case by Sogge and Zelditch [18], [19]. We shall 
use arguments which are similar to those of Shi and Xu [14] and [21] and H. Smith 
(unpublished). 

Specifically, we have the following: 
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Proposition 3.1. If ex is as above and if d = d{x) denotes the distance from x G M to 
dM, 

(3-1) ||(A ^ + d)Vgex\\L°°(M) + \\^x\\l^(m) ^ CX^^\\ex\\Li(^gM)- 

Let us first argue that on the boundary, we have these estimates. Indeed, 

(3.2) A “||IA“eA||Loo(ajvL) < C'aA^ ||eA||Li(aM), 

with here referring to a boundary derivatives. This inequality follows from arguments 
in m and m -[19], since Pex = Xe\ where P is a classical self-adjoint pseudodifferential 
of order one operator whose principal symbol agrees with that of the square root of minus 
the boundary Laplacian. As a result we can use Lemma 5.1.3 in m to write ex = Txex, 
where Tx is an integral operator on the (n — l)-dimensional boundary of M whose kernel 
Kx{x,y) satisfies D'^K = ) for each a, which immediately gives us (|3.2I1 . 

For the next step, we use that by the maximum principle, the bounds in (13.21) for ex 
yield 

n —2 

(3.3) ||eA||L“(M) < C'A^"lleAlliqoM), 

as desired. Thus, we only need to prove the bounds in (|3.1I1 for V^ca. 

As a first step we realize that we can obtain this estimate in the region of M which is of 
distance (5A“^ from the boundary just by using standard Schauder estimates for a given 
<5 > 0. Indeed, since ca is harmonic in M and ()3.2[1 is valid, it follows from Corollary 6.3 
in [8] applied to balls centered at points x £ M or radius r < d{x)l2 that we have 

n —2 

(3.4) ||d VgeA||Lo°({a;eM:dist(x,9M)><5A-i) < C'dA"^" ||eA ||li(9m) ■ 

Here, the constant Cs depends on 6 and {M,g), but not on A. 

To finish the proof of m, it suffices to show that if <5 > 0 is sufficiently small we also 
have the uniform bounds 

■t n — 2 

(3.5) A IIVgeA||L“(MnB(xo,<5A-i)) < C'<5A“2“lleAlliqoM)) xo£dM, 

with B{xo, (5A“^) denoting the geodesic ball of radius (5A“^ about the boundary point xq- 

To prove this we shall use local coordinates and a scaling argument. We shall work in 
such coordinates and scale and normalize ex by replacing it by 

(3.6) uxix) = A ~ex{x/\). 

Similarly, we shall scale the <5A“^ ball so that it becomes a 5 ball B{xo,S) and use the 
“stretched” Laplacian with principal part ^ g^^{x/X)djdk (coming from the “stretched” 
metric gjk{x/X)), which denote by L. It follows from (13.21) that we have the uniform 
bounds 

(3.7) \\D^uxh^^sM)<CJex\\moM), 

where M denotes the stretched version of M in our local coordinates. Additionally, 
the coefficients of our “stretched” Laplacian L belong to a bounded subset of C°° as 
A > 1 and Xg £ dM vary. Also, because of (13.61) we can find a function g}x in our local 
coordinate system which agrees with ux on dM and has bounded C‘^’°‘{B{xo,26) (~l M) 
norm independent of A > 1 and xq £ dM for a given 0 < a < 1. Therefore, if we apply 
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Corollary 8.36 in [S] to u = and / = —Lipx, we conclude that the C^’°‘{B{xo, 6)) 

norm u\ is bonnded uniformly with respect to these parameters if a is fixed. Thus, we 
in particular have the uniform bounds 

\\Dux\\^ 

°°{B(x,S)nM) - 

If we go back to the original local coordinates and recall (EH), we obtain EH, which 
completes the proof of Proposition 13.II 

4. Conclusion. It is now very easy to prove Theorem II. II If we use (j2.6l) and EH, 
we conclude that 

'^lleAllLi(aM) = 2 f I'^exldS < CX^~\\e\\\L^(^gM) [ (A ^ + d(a;)) ^ dS, 

J Z\ Zx 

where, as before, d{x) denotes the distance from x € M to dM. From this, we deduce 
that 

(4.1) {\-^ + d{x))-US. 

Jzx 

Clearly this inequality yields dlH, establishing Theorem ll.il 

Remarks: There is a simple explanation of why the lower bounds (II.5|) are off by a half 
power versus the corresponding best lower bounds dn for the boundaryless case. This 
is because the Dong-type identity in [TB] involved in the left side instead of A, which 
accounts for a relative loss of a full power of A, but, on the other hand, the estimates for 
the gradient here are one half power better due to the fact that the boundary of M is of 
one less dimension, accounting for a relative gain of a half power. 

In some cases one can use (14.11) to get improved lower bounds. For instance if we let 

Zx,k = {x€Zx: d(a:) G [2-^2-'=+l)} 

and if \Zx.k\ < C2~^\Zx\ ioi C <k < log 2 A and if \{x G Zx '■ d{x) < A“^}| < CX~^\Zx\^ 
with C fixed, we then get the lower bound \Zx\ > cA^“' 2 /logA, which is essentially 
optimal when n = 2. The subsets Zx,k of Zx have this property, for instance, for the 
Steklov eigenfunctions r™ sin mO on the disk in (written in polar coordinates). 
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